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Abstract. After briefly recalling some computational aspects of blowing up 
and of representation of resolution data common to a wide range of desin- 
gularization algorithms (in the general case as well as in special cases like 
surfaces or binomial varieties) , we shall proceed to computational applications 
of resolution of singularities in singularity theory and algebraic geometry, also 
touching on relations to algebraic statistics and machine learning. Namely, 
we explain how to compute the intersection form and dual graph of resolution 
for surfaces, how to determine discrepancies, the log-canoncial threshold and 
the topological Zcta-function on the basis of dcsingularization data. Wc shall 
also briefly sec how resolution data comes into play for Bcrnstcin-Sato poly- 
nomials, and wc mention some settings in which dcsingularization algorithms 
can be used for computational experiments. The latter is simply an invitation 
to the readers to think themselves about experiments using existing software, 
whenever it seems suitable for their own work. 



1. Introduction 

This article originated from the notes of an invited talk at the Clay Mathemat- 
ics Institute summer school on " The Resolution of Singular Algcibraic Varieties" in 
Obergurgl, Austria, 2012. As the whole school was devoted to dcsingularization, 
the focus in this particular contribution is on applications and on practical aspcicts. 
A general knowledge of resolution of singularities and different approaches to this 
task is assumed and can be acquired from other parts of this proceedings volume. 
The overall goal of this article is to give readers a first impression of a small choice 
of applications and point them to good sources for further reading on each of the 
subjects. A detailed treatment of each of the topics would fill an article by itself 
and is thus beyond the scope here. 

One focus here is on the practical side. To this end, we first revisit desingularization 
algorithms in section 2 and have a closer look at the representation of resolution 
data: as a consequence of the heavy use of blowing up, the data is distributed over 
a rather large number of charts which need to be glued appropriately. Glueing, 
however, only describes the theoretical side of the process; from the practical point 
of view, it is closer to an identification of common points in charts. 
In section 3 we focus on applications needing different amounts of resolution data. 
Using an abstract resolution of singularities only, the computation of the inter- 
section form and dual graph of the resolution for surface singularities requires 
the smallest amount of data. For determining discrepancies and the log-canonical 
threshold, we already need an embedded resolution which is also required for the 
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third application, the computation of the topological zeta function. 
In the rather short last section, we only sketch two settings in which one might 
want to use algorithmic desingularization as an experimental tool: the roots of 
the Bernstein-Sato polynomial and resolution experiments in positive characteris- 
tic. This last part is not intended to provide actual research projects. It is only 
intended to help develop a feeling for settings in which experiments can be helpful. 
All parts of the article are illustrated by the same simple example which is desin- 
gularized using a variant of Villamayor's algorithm available in Singular. Based 
on this resolution data, all further applications are also accompanied by the cor- 
responding SiNGULAR-code. The Singular code is not explained in detail, but 
hints and explanations on the appearing commands and their output are provided 
as comments in the examples. To further familiarize with the use of Singular in 
this context, we recommend that the readers try out the given session themselves 
and use the built-in manual of Singular to obtain further information on the com- 
mands (e.g.: help resolve; returns the help page of the command resolve). 
I would like to thank the organizers of the summer school for the invitation. In- 
sights from conversations with many colleagues have contributed to the content 
of these notes. In particular, I would like to thank Herwig Hauser, Gerhard Pfis- 
ter, Ignacio Luengo, Alejandro Melle, Frank-Olaf Schreyer, Wolfram Decker, Hans 
Schonemann, Duco van Straten, Nobuki Takayama, Shaowei Lin, Frank Kiraly, 
Bernd Sturmfels, Zach Teitler, Nero Budur, Rocio Blanco and Santiago Encinas 
for comments, of which some led to applications explained here and some others 
helped seeing the applications in broader context. For reading earlier versions of 
this article and many helpful questions on the subject of this article, I am indepted 
to Frithjof Schulze and Bas Heijne. 

2. Desingularization from the computational side 

Before turning toward applications of resolution of singularities, we need to re- 
view certain aspects of algorithmic desingularization to understand the way in which 
the computed resolution data is represented. Although there are various settings 
in which different resolution algorithms have been created, we may discern three 
main approaches suitable for the purposes of this article: the algorithms based on 
Hironaka's proof in characteristic zero in any dimension (see e.g. [1], [S], [13]), the 
algorithms for binomial and toric ideals (see e.g. [15], [7], [6]) and the algorithms 
for 2-dimensional varieties and schemes (such as [5] - based on [T^] - or [24i)Q The 
algorithms of the first kind of approach involve embedded desingularization, i.e. 
they blow up a smooth ambient space and consider the strict transform of the vari- 
ety and the exceptional divisors inside the new ambient space. The algorithms for 
2-dimensional varieties on the other hand, do not consider the embedded situation, 
but blow up the variety itself and consider exceptional divisors inside the blown up 
variety. 

It would be beyond the scope of this article to cover all these algorithms in 
depth, but they all have certain ingredients in common. In the first two situations. 
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a desingularization is achieved by finite sequences of blow-ups at suitable non- 
singular centers. The differences between these algorithms then lie in the choice of 
center, which is the key step of each of these, but does not affect the structure of 
the practical representation of resolution data. For the third class of algorithms, 
blow-ups are not the only tool, but are combined with other tools, in particular 
normalization steps. However, the exceptional divisors to be studied arise from 
blow-ups and additionally only require proper tracing through the normalization 
steps if necessary. Therefore the technique to focus on in this context is blowing 
up; more precisely blowing up at non-singular centers. 

2.1. Blowing up — the computational side. Let us briefly recall the definition 
of blowing up, as it can be found in any textbook on algebraic geometry (e.g. [IS]), 
before explaining its computational side: 

Definition 1. Let X be a scheme and Z G X a suhscheme corresponding to a 
coherent ideal sheaf T. The blowing up of X with center Z is 

^ -.X := Proji^I^) X. 

d>0 

Let Y^X be a closed suhscheme and tti : Y — > Y the blow up ofY along i^^IOy ■ 
Then the following diagram commutes 

Y X 

Y X 

Y is called the strict transform ofY, 7t*{Y) the total transform of Y . 

To make X accessible to explicit computations of examples in computer algebra 
systems, it should best be described as the zero set of an ideal in a suitable ring. 
We shall assume now for simplicity of presentation that X is affine because schemes 
are usually represented in computer algebra systems by means of affine covers. So 
we are dealing with the following situation: J — (/i, . . . , f„i) C A is the vanishing 
ideal of the center Z <Z X = Spec{A) and the task is to compute 

d>0 

To this end, we consider the canonical graded A-algebra homomorphism 
$:A[2/i,...,y„]— >0jVncA[i] 

n>0 

defined by 3>(yi) = tfi. The desired object ®(j>o is then isomorphic to 

A[yi,...,y^]/ker{^) 
or from the more geometric point of view V{ker{^)) C Spec{A) x P'""-'^. 

The computation of the kernel in the above considerations is a standard basis 
computation and further such computations arise during the calculation of suit- 
able centers in the different algorithms. Moreover, each blowing up introduces a 
further set of new variables as we have just seen and desingularization is hardly 
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ever achieved with just one or two blow-ups - usually we are seeing long sequences 
thereof. The performance of standard bases based algorithms, on the other hand, 
is very sensitive to the number of variables as its complexity is doubly exponential 
in this number. Therefore it is vital from the practical point of view to pass from 
the P™~i to TO afhne charts and pursue the resolution process further in each of 
the chart^ consistency of the choice of centers does not pose a problem at this 
stage as this follows from the underlying desingularization algorithm. Although 
this creates an often very large tree of charts with the final charts being the leaves 
and although it postpones a certain part of the work, this approach has a further 
important advantage: it allows parallel computation by treating several charts on 
different processors or cores at the same time. 

As a sideremark, we also want to mention the computation of the transforms, be- 
cause they appeared in the definition cited above; for simplicity of notation we only 
state the afhne case. Without any computational effort, we obtain the exceptional 
divisor as /(-ff) = JO^ and the total transform 7r*(/) = lO^ for a subvariety 
ViV) hi our afhne chart. The strict transform is then obtained by a saturation, i.e. 
an iteration of ideal quotients until it stabilizes: ^^[1) ~ (^*('^^ ■ for the 

weak transform, the iteration stops before it stabilizes, namely at the point where 
multiplying with the ideal of the exceptional divisor gives back the result of the 
previous iteration step for the last time. These ideal quotient computations are 
again based on standard bases. 

2.2. Identification of Points in Different Charts. As we have just seen, it is 
more useful for the overall performance to pass to affine charts after each blow-up, 
even though this leads to an often rather large tree of charts. As a consequence, 
it is not possible to directly work with the result without any preparation steps, 
namely identifying points which are present in more than one chart - a practical 
step equivalent to the glueing of the charts. This is of particular interest for the 
identification of exceptional divisors which are present in more than one chart. 

For the identification of points in different charts, we need to pass through the 
tree of charts - from one final chart all the way back to the last common ancestor 
of the two charts and then forward to the other final chart. As blowing up is an 
isomorphism away from the center, this step does not pose any problems as long 
as we do not need to identify points lying on an exceptional divisor which was not 
yet created in the last common ancestor chart. In the latter case, however, we do 
not have a direct means of keeping track of points originating from the same point 
of the center. The way out of this dilemma is a representation of points on the 
exceptional divisor as the intersection of the exceptional divisor with an auxilliary 
variety not contained in the exceptional divisor. More formally, the following simple 
observation from commutative algebra can be used: 

Remark 2. Let I C K[xi, . . . , x^] be a prime ideal, J C K[ai, . . . , Xn] another 
ideal such that I + J is equidimensional and ht{I) — ht{I + J) ~ r for some integer 
< r < n. Then there exist polynomials pi, . . . ,pr G I + J and a polynomial 



In practice, it is very useful to discard all charts not appearing in any other chart and not 
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/ e K[xi, . . . , Xn] such that 



VTTJ=^iI+{pi,...,Pr)):f. 



At first glance, this seems non-constructive, but it turns out to be applicable 
in a very convenient way: In our situation, / describes the intersection of ex- 
ceptional divisors containing the point (or subvariety V{J)) in question. As any 
sufficiently general set of polynomials pi, . . . ,pr S J\ {I D J) with correct height 
of / + (pi, . . . ,pr) will serve our purpose and so will any / that excludes all extra 
components of / + {pi, . . . ,Pr), we have sufficient freedom of choice of pi, . . . ,pr, f 
such that none of these is contained in any further exceptional divisor whose mo- 
ment of birth has to be crossed during the blowing down process. 

Given this means of identification of points, we can now also identify the excep- 
tional divisors or, more precisely, the centers leading to the respective exceptional 
divisors. To avoid unnecessary comparisons between centers in different charts, we 
can a priori rule out all comparisons involving centers lying in different exceptional 
divisors. If the desingularization is controlled by an invariant as in [5], [5] or [7], 
we can also avoid comparisons with different values of the controlling invariant, 
because these cannot give birth to the same exceptional divisor either. 

Example 3. To illustrate the explanations given so far and to provide a practical 
example to be used for all further applications, we now consider an isolated surface 
singularity of type at the origin. We shall illustrate this example using the 
computer algebra system SINGULAR ( (llj )■ 

> // load the appropriate libraries for resolution of singularities and 

> // applications thereof 

> LIB"resolve.lib" ; 

> LIB"reszeta.lib" ; 

> LIB"resgraph.lib" ; 

> // define the singularity 

> ring R=0, (x,y,z) ,dp; 

> ideal I=x5+y2+z2; // sm. A4 surface singularity 

> // compute a resolution of the singularity (Villamayor-approach) 

> list L=resolve(I) ; 

> size(L[l]); // final charts 
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> size (L [2] ) ; 
11 



// all charts 



> def r9=L[2] [9] ; 



// go to chart 9 



> setring r9; 

> showBD(BD) ; 



// show data in chart 9 



==== Ambient Space: 
_[1]=0 



==== Ideal of Variety: 
_[l]=x(2)-2+y(0)+l 
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==== Exceptional Divisors: 



[1] : 

_[!]=! 
[2] : 

_[l]=y(0) 
[3] : 

_[!]=! 
[4] : 

_[l]=x(l) 

==== Images of variables of original ring: 
_[l]=x(l)-2*y(0) 
_[2]=x(l)-5*x(2)*y(0)-2 
_[3]=x(l)-5*y(0)-2 

> setring R; // go back to old ring 

This yields a total number of 11 charts of which 6 are final charts. All blow-ups 
have zero- dimensional centers except the two line blow-ups leading from chart 6 to 
charts 8 and 9 and from chart 7 to charts 10 and 11. As it would not be very 
useful to reproduce all data of this resolution here, we show the total tree of charts 
as illustration^^ and give the content of chart 9 as an example: 



y I — !■ xlx2yo 
z I — > x\yl 

For example, the exceptional divisor originating from the last blow-up leading to 
chart 9 needs to be compared to the exceptional divisor originating from the last blow- 
up leading to chart 10. To this end, one needs to consider the centers computed in 
charts 6 and 7 which turn out to be the intersection of the exceptional divisors labeled 
2 and 3, if one considers the output of the resolution process in detai^ Therefore 
the last exceptional divisors in charts 9 and 10 coincide. This identification is 
implemented in SINGULAR and can be used in the following way: 

> // identify the exceptional divisors 

> list coll=collectDiv(L) ; 

> coll[l] ; 

0,0,0,0, // no exc . div. in chart 1 

1,0,0,0, // first exc. div. is first in chart 2 

1,0,0,0, // 

1,2,0,0, // data too hard to read, better 
0,2,0,0, // use command below to create figure 1 




V{xi) from last blow-up 



images of variables of original ring: 



x I — > xjyo 



Weak and strict transform coincide for hypcrsurfaccs. 
^We encourage the reader to verify this by typing the above sequence of commands into SIN- 
GULAR and then exploring the data in the different charts. 
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chai1 2 




chaj1 3 


E:l 




E:l 




chail 4 




c\vM^ 5 


E:l,2 




E:2 




chart 6 

E:2.3 




chajl 7 

E:2,3 










chati 9 
E:2,4 




chart 10 
E:3,4 



Figure 1 . Tree of charts of an embedded desingularization of an 
A4 surface singularity. The numbers given in the second hne in 
each chart arc the labels of the exceptional divisors visible in the 
respective chart. The numbers stated as d = or d = 1 state 
the dimension of the center of the corresponding blow-up. Charts 
providing only data which is also present in other charts are not 
shown. 



0,2,3,0, 
0,2,3,0, 
0,0,3,4, 
0,2,0,4, 
0,0,3,4, 
0,2,0,4 



> //present the tree of charts as shown in figure 1 

> ResTree(L,coll[l]) ; 

3. Applications of Resolution of Singularities 

The applications we present in this section originate from different subfields of 
mathematics ranging from algebraic geometry to singularity theory and D-modules. 

For each application we shall revisit our example from the previous section and also 
show how to perform the corresponding computation using SINGULAR. 

3.1. Intersection Form and Dual Graph of Resolution. Given a resolution 
of an isolated surface singularity, we want to compute the intersection matrix of 
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the excpetional divisors. This task does not require an embedded resolution of 
singularities, only an abstract one. Given such a desingularization, it can then be 
split up into 3 different subtasks: 

(1) computation of the intersections Ei.Ej for exceptional curves Ei ^ Ej 

(2) computation of the self-intersection numbers Ef for the exceptional curves 

E^ 

(3) representation of the result as the dual graph of the resolution 

If the given resolution, is not an abstract one, but an embedded one - like the 
result of Villamayor's algorithm - we need to add a preliminary step 

(0) determine an abstract resolution from an embedded one. 
Although the definition of intersection numbers of divisors on surfaces can be 
found in many textbooks on algebraic geometry (e.g. in [TB], V.l), we give a brief 
summary of the used properties for readers' convenience: 

Definition 4. Let Di, D2 be divisors in general positiot^ on a non-singular surface 
X. Then the intersection number of Di and D2 is defined as 

D,.D2 := (^i-^2)x 

xeDinD2 

where {Di.D2)x denotes the intersection multiplicity of Di and D2 at x. 

Lemma 5. For any divisors Di and D2 on a non-singular surface X , there exist 
divisors D[ and D'2, linearly equivalent to Di and D2 respectively, such that D[ 
and D'2 are in general position. 

Lemma 6. Intersection numbers have the following basic properties: 

(a) For any divisors C and D: CD = D.C . 

(b) For any divisors C, Di and D2: C.{Di + D2) = C.Di + C.D2 

(c) For any divisors C, Di and D2, such that Di and D2 are linearly equivalent: 
C.Di ^ C.D2. 

At this point, we know what we want to compute, but we have to take care 
of another practical problem before proceeding to the actual computation which is 
then a straight forward calculation of the intersection numbers of exceptional curves 
Ei ^ Ej . The practical problem is that computations in a computer algebra system 
usually take place in polynomial rings over the rationals or algebraic extensions 
thereof. So we easily achieve a decomposition of the exceptional divisor into Q- 
irreducible components, but we need to consider C-irreducible components to obtain 
the intersection matrix we expect from the theoretical point of view. To this end, 
passing to suitable extensions of the ground field may be necessary - be it explicitly 
by introducing a new variable and a minimal polynomial (slowing down subsequent 
computations) or implicitly by taking into account the number of components over 
C for each Q-component. 

Example 7. Revisiting our example of a desingularization of an A4 surface sin- 
gularity, we first need to pass to an abstract resolution. 

^This is achieved by canceling all trailing blow-ups in our tree of charts which are unnecessary 
for the non-embedded case. Then the intersection of the remaining excpetional divisors with the 
strict transform yields the exceptional locus of the non-embedded resolution. 

^Di and D2 are in general position, if the intersection Supp{Di) n Supp{D2) is either empty 
or a finite set of points. 
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\\ compute part of the tree of charts relevEint for abstract resolution 

> abstractR(L) [1] ; 

0,1,0,1,1,0,0,0,0,0,0 //final charts are 2,4,5 

> abstractR(L) [2] ; 

0,0,0,0,0,1,1,1,1,1,1 //charts 6 and higher are irrelevant 

//for non-embedded case 

So we only see 2 exceptional divisors in the final charts, the ones labeled 1 and 2 ( cf. 
figure^. But looking at the charts in more detail, we would see e.g. in chart 4 that 
the first divisor is given by V{y2,yQ + 1) and the second one by V{x2,yQ + 1). So 
each of these has two C-irreducible components. Considering these 'C- components, 
we can obtain the following intersection data (seeing two of the intersections directly 
in chart 4 and the remaining one in chart 5): 

/ * 1 \ 

0*01 

10*1 
\ 1 1 * / 

Hence only the self-intersection numbers - marked as * in the matrix - are still 
missing. 

For the self-intersection numbers of the exceptional curves, we need to make use 
of another property of divisors in the context of desingularization; 

Lemma 8. Let n : X — > X be a resolution of singularities of a surface X . Let 
Di be a divisor on X all of whose components are exceptional curves of tt and let 
D2 be a divisor on X , then 

Denoting hy Ei, . . . , Eg the C-irreducible exceptional curves, we can hence con- 
sider a linear form h : X — > C passing through the only singular point of X and 
the divisor D defined by it. We then know 

s 

^^*{D)=Y,c^E, + H 

4 = 1 

where H denotes the strict transform of D and the Ci are suitable integers. From 
the lemma we additionally know that 

= TT*{D).Ej = c,E,.Ej + H.Ej VI < j < s 

1=1 

where all intersection numbers are known or directly computable in each of the 
equations except the self-intersection numbers Ej.Ej which we can compute in this 
way. For the dual graph of the resolution, each divisor is represented by a vertex 
(those with self-intersection -2 are unlabeled, the other ones labeled by their self- 
intersection number), each intersection is represented by an edge linking the two 
vertices corresponding to the intersecting exceptional curves. 

Example 9. The computation of the intersection form is implemented in SINGU- 
LAR and can be used as follows: 

> 1 1 intersection matrix of exceptional curves 

> // (no previous abstractR is needed, this is done automatically 
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> // in the procedure intersectionDiv) 

> list iD=iiitersectionDiv(L) ; 

> iD[l]; 
-2,0,1,0, 
0,-2,0,1, 
1,0,-2,1, 
0,1,1,-2 

> // draw the dual graph of the resolution 
InterDiv(iD[l] ) ; 

This yields the expected intersection matrix (as entry iD[lJ of the result) 

( -2 1 Q\ 
0-201 
10-21 
\ 1 1 -2 / 

and the dual graph of the resolution which is just the Dynkin diagram of the A4 
singularity. 

3.2. Discrepancies and Log-canonical Threshold. In contrast to the last task, 
which only required an abstract resolution of the given surface singularity, the task 
of computing (log-)discrepancies and the log-canonical threshold requires embedded 
desingularization (or principalization of ideals). This is provided by Villamayor's 
algorithm. As before we first recall the definitions and some basic properties (see e.g. 
[25] for a direct and accessible introduction to the topic) . To keep the exposition as 
short as possible and the considerations directly accesible to explicit computation, 
we restrict our treatment here to the case of a singular affine variety. 

Definition 10. Let f G C[a;i, .... Xn] he a non-zero polynomial defining a hyper- 
surface V and let TT : X — > C" be an embedded resolution of V . Denote by Ei, 
i£ I, the irreducible components of the divisor iT~'^{f~^{0)). Let N(Ej) denote the 
multiplicity of Ej, j € J in the divisor of f ot: and let v{Ej) — 1 he the multiplicity 
of Ej in the divisor Kx/c^ — TT*{dxi A ... A rfa;„). Then the log-discrepancies of the 
pair (C",y) w.r.t. Ej, j € J, are 

a{Ej;U\ V) := v{Ej) - N{Ej). 

The minimal log discrepancy of the pair (C", V) along a closed subset W C C" is 
the minimum over the log-discrepancies for all Ej with TT{Ej) C W , i.e. originating 
from (sequences of) blow-ups with centers in W . The log-canonical threshold of the 
pair (C", V) is defined as 



lct{'C-,V)^inf,^., 



HE,) 
N{E,y 



Remark 11. The above definition of log-discrepancies and log-canonical threshold 
holds in a far broader context. Allowing more general pairs (Y, V) it is also the basis 
for calling a resolution of singularities log- canonical, if the minimal log discrepancy 
of the pair along all of Y is non-negative, and log-terminal, if it is positive. 

As we already achieved an identification of exceptional divisors in a previous 
section, the only computational task here is the computation of the multiplicities 
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N{Ei) and The fact that we might be deahng with Q-irreducible, but C- 

reducible Ei does not pose any problem here, because we can easily check that the 
respective multiplicities coincide for all C-components of the same Q-component. 
To compute these multiplicities from the resolution data in the final charts (i.e. 
without moving through the tree of charts), we can determine N{Ei) by finding the 
highest exponent j such that I{EjY : J is still the whole ring where J denotes the 
ideal of the total transform of the original variety. To compute the i'{Ei) we can use 
the same approach, but taking into account the appropriate Jacobian determinant. 

Example 12. We now simply continue our SINGULAR session on the basis of the 
data already computed in the previous examples 

>// identify exceptional divisors (embedded case) 

> list iden=prepEmbDiv(L) ; 

>// multiplicities N(Ei) 

> intvec cN=computeN(L,iden) ; 

> cN; // last integer is strict transform 
2,4,5,10,1 

>// multiplicities v(Ei) 

> intvec cV=computeV(L,iden) ; 

> cV; // last integer is strict transform 
3,5,7,12,1 

>// log-discrepancies 

> discrepancy (L) ; 
0,0,1,1 

>// compute log-CEinonical threshold 
>// as an example of a loop in Singular 

> number lct=number(cV[l] )/number(cN[l] ) ; 

> number Icttemp; 

> for (int i=l; i < size(cV); i++) 

> { 

> lcttemp=number(cV[i] ) /number (cN[i] ) ; 

> if (Icttemp < let) 

> { 

> lct=lcttemp; 

> } 

> } 

> let; 
6/5 

The log-canonical threshold, in particular, is a very important invariant which 
appears in many different contexts, ranging from a rather direct study of proper- 
ties of pairs to the study of multiplier ideals (cf. [32]), to motivic integration or to 
Tian's a-invariant which provides a criterion for the existence of Kahler-Einstein 
metrics (cf. [25] )■ 
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A real analogue to the log-canonical threshold, the real log-canonical threshold 
appears when applying resolution of singularities in the real setting [55]. In alge- 
braic statistics, more precisely in model selection in Bayesian statistics, desingular- 
ization plays an important role in understanding singular models by monomializing 
the so-called KuUback-Leibler function at the true distribution. In this context the 
real log-canonical threshold is then used to study the asymptotics of the likelihood 
integral. [23j. It also appears in singular learning theory as the learning coefficient. 

3.3. Topological Zeta Function. Building upon the multiplicities of N{Ei) and 
v{Ei) — 1 of exceptional curves appearing in / o tt and Kx/C" — ^* (dxi A ... A dxn) 
which already appeared in the previous section, we can now define and compute 
the topological Zeta-function. As before, we first recall the definitions and prop- 
erties (see e.g. p5]) and then continue with the computational aspects - listing 
the corresponding Singular commands in the continued example of an A4 surface 
singularity. 



Definition 13. Let f G C[a:i, . . . , Xn] be a non-zero polynomial defining a hyper- 
surface V and let TT : X — > C" be an embedded resolution of V . Denote by Ei, 
i G I, the irreducible components of the divisor 7r^^(/^^(0)) . To fix notation, we 
define for each subset J C I 

Ej := rijejEj and E*j := Ej \ Uj^jEj^ij} 
and denote for each j (z I the multiplicity of Ej in the divisor of f o t: by N{Ej). 
We further set ^{Ej) — 1 to be the multiplicity of Ej in the divisor Kx/c^ — 
'K*(dxi A ... A dxn)- In this notation the topological Zeta-function of f is 

z'itlU,^) E x{E:j)\{HE,)+N{Ej)s)-^ G^{s). 

JC/a.th. j£j 
d\N{Ej)VjeJ 

Intersecting the Ej with the preimage of zero in the above formula leads to the local 
topological Zeta-function 

JC-fB.th. ,-gJ 

(The local and global topological Zeta-function are independent of the choice of 
embedded resolution of singularities of V .) 

Here, it is again important to observe that in the above context the irreducible 
components are taken over C, while practical calculations usually take place over 
Q and further passing to components taken over C is rather expensive. The fol- 
lowing lemma shows that considering Q-irreducible components already allows the 
computation of the ^-function: 

Lemma 14. Let Di ,IgL, be the Q-irreducible components of the divisor Tr~^{f^^{0)). 
For each subset J d L define Dj and Dj as above. Then the topological (global 
and local) Zeta-function can be computed by the above formulae using the Dj and 
Dj instead of the Ej and Ej . 

As we already identified the exceptional divisors and computed the multiplicities 
N[Ei) and i^{Ei), the only computational data missing is the Euler characteristic 
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of the exceptional components in the final charts. If the intersection of excep- 
tional divisors is zero-dimensional, this is just a matter of counting points using 
the identification of points in different charts. For 1-dimensional intersections the 
Euler characteristic can be computed using the geometric genus of the curve (using 
x(C) = 2 — 2g(C)). Starting from dimension two on, this becomes more subtle. 

Example 15. In our example which we have been treating throughout this article, 
we are dealing with a surface in a three-dimensional ambient space. So the only 
further Euler characteristics which need to be determined are those of the exceptional 
divisors themselves. At the moment of birth of an exceptional divisor, it will either 
be a with Euler characteristic 3 or a x C (for a one- dimensional center C ) 
leading to Euler characteristic 4 — 'ig{C). Under subsequent blow-ups the tracking 
of the changes to the Euler characteristic is then no difficult task. 
II compute the topological zeta-f unction for our isolated 
// surface singularity (global Eind local zeta-f unction coincide) 

> zetaDL(L,l); // global zeta-f unction 
[1] : 

(s+6)/(5s2+lls+6) 

> zetaDL(L, 1 , "local") ; // local zeta-f unction 
Local Case: Assuming that no (!) charts were dropped 
during calculation of the resolution (option "A") 

[1] : 

(s+6)/(5s2+lls+6) 

// zetaDL also computes the characteristic polynomial 

// of the monodromy, if additional parameter "A" is given 

> zetaDL(L,l,"A") ; 
Computing global zeta function 
[1] : 

(s+6)/(5s2+lls+6) 
[2] : 

(s4+s3+s2+s+l) 

4. DESINGULARIZATION IN EXPERIMENTS 

The previous section showed some examples in which desingularization was a 
crucial step in the calculation of certain invariants and was hence used as a theo- 
retical and practical tool. We now turn our interest to a different kind of settings: 
experiments on open questions which involve desingularization. Here we only sketch 
two such topics and the way one could experiment in the respective setting^ 

4.1. Bernstein-Sato polynomials. In the early 1970s J.Bernstein [1] and M.Sato 

[2 7) independently defined an object in the theory of D-modulcs, which is nowadays 
called the Bernstein-Sato polynomial. The roots of such polynomials have a close, 
but still somewhat mysterious relation to the multiplicities of exceptional divisors 



Neither of the two topics should be seen as a suggestion for a short term research project! Both 
questions, however, might gain new insights from someone playing around with such experiments 
just for a short while and stumbling into examples which open up knew perspectives, insight or 
conjectures. 
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in a related desingularization. For briefly recalling the definition of Bcrnstcin-Sato 
polynomials, we shall follow the article of Kashiwara [21 , which also introduces 
this relation. After that we sketch what computer algebra tools are available in 
Singular for experiments on this topic. 

Definition 16. Let f be an analytic function defined on some complex manifold 
X of dimension n and let V be the sheaf of differential operators of finite order on 
X. The polynomials in an additional variable s satisfying 

form an ideal. A generator of this ideal is called the Bernstein- Sato polynomial and 
denoted by bf{s). 

Kashiwara then proves the rationality of the roots of the Bernstein-Sato polyno- 
mial by using Hironaka's desingularization theorem in the following way: 

Theorem 17 ([21 ). Let f be as above and consider a blow-up F : X' — X and a 
new function f — foF. Then bf{s) is a divisor o/J^^g bf'{s + k) for a sufficiently 
large N . 

More precisely, a principalization of the ideal generated by / leads to a poly- 
nomial /' = Y[i=i '^^r with a local system of coordinates ti, . . . ,tm- For /' the 
Bernstein-Sato polynomial is known to be 

m Ti 

bf'{s)^l[l[{r,s-^k). 

i=i k=i 

Therefore the roots of the Bernstein-Sato polynomial are negative rational numbers 
of the form — ^ with exceptional mulitplicities and 1 < k < ri. Yet it is not 
clear at all whether there are any rules or patterns which exceptional multiplicities 
ri actually appear as denominators of roots of the Bernstein-Sato polynomial. (As 
a sideremark: If this were known, then such knowledge could be used to speed 
up computations of Bernstein-Sato polynomials whenever a desingularization is 
known.) 

Today there are implemented algorithms for computing the Bernstein-Sato poly- 
nomial of a given / available: by Ucha and Castro- Jiminez [29], Andres, Levan- 
dovskyy and Martin- Morales [T] and by Berkesh and Leykin [3]. All of these algo- 
rithms do not use desingularization techniques, but rather rely on Grobner-Bases 
and annihilator computations. The computed objects, however, are also used (the- 
oretically and in examples) in the context of multiplier ideals and thus have relation 
to invariants like the log-canonical threshold. 

Given algorithmic ways to independently compute the Bcrnstein-Sato polynomial 
and the exceptional multiplicities, one could now revisit the exploration of the 
interplay between these and try to spot patterns to get a better understanding. 

4.2. Positive Characteristic. Desingularization in positive characteristic is one 
of the long standing, central open problems in algebraic geometry. In dimensions 
up to three there is a positive answer (see e.g. [5], [IH]), but in the general case 
there are several different approaches (e.g. [20], [30], [IH]) each of which has run 



DESINGULARIZATION IN COMPUTATIONAL APPLICATIONS AND EXPERIMENTS 15 



into obstacles which are currently not resolved. 

About a decade ago, Hauser started studying the reasons why Hironaka's ap- 
proach of characteristic zero fails in positive characteristic [17j . Among other find- 
ings, he singled out two central points which break down: 

(1) failure of maximal contact: 

Hypersurfaces of maximal contact are central to the descent in ambient di- 
mension which in turn is the key to finding the correct centers for blowing 
up. In positive characteristic, it is well known that hypersurfaces of max- 
imal contact need not exist; allowing hypersurfaces satisfying only slightly 
weaker conditions is one of the central steps in the approach of Hauser and 
Wagner for dimension 2 |18| . (For higher dimensions this definition requires 
a little bit more care [H].) 

(2) increase of order of coefficient ideal: 

As the improvement of the singularities is measured by the decrease of 
the order and of orders for further auxilliary ideals constructed by means 
of descent in ambient dimension, it is crucial for the proof of termination 
of resolution that these orders cannot increase under blowing up. Unfor- 
tunately this does no longer hold for the orders of the auxilliary ideals 
in positive characteristic as has again been known since the 1970s. Hauser 
characterized the structure of polynomials which can exhibit such behaviour 
in [17]. 

Although problems of desingularization in positive characteristic are known, this 
knowledge seems to be not yet broad enough to provide sufficient feedback for suit- 
able modification of one of the approaches to overcome the respective obstacles. 
Experiments could prove to be helpful to open up a new point of view. In par- 
ticular, the approach of Hauser and Wagner for surfaces is sufficiently close to the 
characteristic zero approach of Hironaka (and hence to algorithmic approaches like 
the one of Villamayor) to allow modification of an existing implementation to pro- 
vide a tool for a structured search for examples with special properties also in higher 
dimensions. This has e.g. been pursued in [T3]. 
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